Tuning the conductance of molecular junctions: transparent versus tunneling regimes 
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We present a theoretical study of the transport characteristics of molecular junctions, where first-row diatomic 
molecules are attached to (001) gold and platinum electrodes. We find that the conductance of all of these 
junctions is of the order of the conductance quantum unit Go, spelling out that they belong to the transparent 
regime. We further find that the transmission coefficients show wide plateaus as a function of the energy, instead 
of the usual sharp resonances that signal the molecular levels in the tunneling regime. We use Caroli's model 
to show that this is a rather generic property of the transparent regime of a junction, which is driven by a strong 
effective coupling between the delocalized molecular levels and the conduction channels at the electrodes. We 
analyse the transmission coefficients and chemical bonding of gold/Benzene and gold/Benzene-dithiolate (BDT) 
junctions to understand why the later show large resistances, while the former are highly conductive. 

PACS numbers: 73.63.Rt,73.40.-c,73.63.-b 



INTRODUCTION 



The field of molecular electronics was arisen by the earl 
realization that organic molecules could act as rectifiers jl 
when attached to conducting electrodes to form tunnel junc- 
tions. Many experiments with a large variety of organic 
molecules have been performed!!! 0], typically finding val- 
ues of the conductance G several orders of magnitude smaller 
than Go (Go = 2e 2 /h is the conductance quantum) and a 
large variability, which hinder the reproducibility of the ex- 
periments. Molecular junctions can be understood in terms 
of resonant tunneling modelsj^l, where the conduction is car- 
ried through the Highest Occupied and Lowest Unoccupied 
Molecular Orbitals (HOMO and LUMO, respectively). These 
are revealed as sharp resonances in either the Densities of 
States (DOS) of the molecule, or the Transmission coefficients 
T(E) of the junction, and are usually located 1 or 2 eV above 
or below the Fermi level of the molecule, respectively. Con- 
ductance values of the order of Go can only be achieved by 
pinning one of those resonances to the Fermi level of the elec- 
trodes. Otherwise, the conductance is very low. 

The invention of the Scanning Tunneling Microscopel|5|l al- 
lowed the fabrication of stable atomic point contacts J(| 0, 0, 
I9I [Toll . These junctions were found to be highly transparent 
in many cases, and to show values of the conductance of the 
order of Gn, which confirmed early theoretical predictions on 
the matter ! 1 ll Il2ll . Theoretical analyses of these junctions 



Importantly, more recent developments using MCBJ tech- 
niques demonstrated that high values of the conductance (~ 
Go) are not restricted to atomic constrictions, but could also 
be obtained even when platinum or palladium electrodes are 



found that their transmission coefficients T(E) show wide 
plateaus as a function of the energy E of the incoming elec- 
trons, with heights of order one. The high transparency of 
these junctions is due to the good matching between the con- 
duction channels at the electrodes and those at the molecule. 
The (contact) resistance of the junction is different from zero 
because of the different number of channels at the electrodes 
and the junction which leads to a recombination of the former 
to match the later fl [Till. 



bridged by hydrogen molecules ! 1 31 1 1411. T hese results were 
reproduced by theoretical simulations rtl5l which also 
showed how the antibonding level of the hydrogen molecule 
hybridized strongly with the conduction channels of the elec- 
trodes, providing a junction with a single conduction channel. 
The transparency of this channel was manifested in the trans- 
mission coefficients T(E), that had a wide plateau of height 
one. Furthermore, very recent experimental work confirms 
that junctions comprising platinum electrodes and either sim- 
ple benzene! 1711, or a number of small molecules] 18] show 
conductance values of the order of Go . 

As stated above, junctions in the tunneling regime can shed 
conductance values of order Go provided that a molecular 
level is exactly pinned to the Fermi energy. But it is hard 
to believe that all of the junctions discussed in those recent 
experiments! T7I 18 1 display this pinning mechanism. Instead, 
they clearly indicate that highly transparent molecular junc- 
tions can be fabricated with relative ease. They also indirectly 
hint that thiol-capping necessarily leads to junctions in the 
tunneling regime. 

We have performed a number of transport simulations of 
molecular junctions where first-row diatomic molecules are 
sandwiched by semi-infinite (001 ) g old and platinum elec- 
trodes, with our code SMEAG0L !19I1 . Our simulations con- 
firm that this type of junctions is highly conductive. Indeed, 
the transmission coefficients do not show resonant behavior 
but wide plateaus of height of order one instead. To sus- 
tain theoretically these simulations, we argue that the con- 
ductive behavior of a junction is determined by two factors. 
First, by the conjugation nature of the molecule, e.g.: whether 
the HOMO and LUMO levels are delocalized throughout the 
whole molecule, or not. Second, by the strength of the chem- 
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FIG. 1: Transmission coefficients at zero voltage as a function of 
energy T(E), for a number of first-row diatomic molecules and CO 
attached to (001) gold electrodes. The zero of energies corresponds 
to the position of the Fermi energy of the electrodes at zero voltage. 
The distance between the last gold atoms at the electrodes and the 
molecules is set to the equilibrium distance for the gold-i?2 junction 
(1.5 A). 
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FIG. 2: Transmission coefficients at zero voltage as a function of 
energy for a number of first-row diatomic molecules and CO attached 
to (001) platinum electrodes. The geometry of the junctions has been 
relaxed in this case. 



SIMULATIONS OF FIRST-ROW DIATOMIC MOLECULES 
CONTACTED TO (001) GOLD AND PLATINUM 
ELECTRODES 



ical bond between the conduction channels at the electrodes 
and these delocalized HOMO or LUMO orbitals. A junction 
will be highly conductive provided that its molecule is conju- 
gated and the chemical bond referred above is strong. 

Conjugated molecules are reasonably well approximated by 
Caroli's model I20I1 . We use this model below to show that 
when the chemical bond between the delocalized HOMO or 
LUMO orbitals and the conduction channels is strong, the 
transmission coefficients show wide plateaus whose height is 
of the order of Go- We argue that these plateaus are robust 
against changes in the energy of the HOMO/LUMO levels, so 
that there is no need to fine-tune their position to the Fermi 
level of the electrodes in order to achieve large conductance 
values. On the contrary, when the conduction channels at the 
electrodes do not bind chemically with neither the HOMO nor 
the LUMO levels then the junction is in the tunneling regime 
where it shows resonant behavior. 

Diatomic or triatomic molecules are sufficiently small to 
show conjugation. Likewise, Benzene-based molecules are 
archetypical conjugated molecules, where the conjugation is 
driven by internal 7r-bonding. It is therefore important to un- 
derstand unequivocally why BDT junctions show strong tun- 
neling behavior while on the contrary the simpler Benzene 
junctions are highly transparent. We perform below a study 
of the transmission coefficients of these junctions. Our study 
indicates that direct Au-S bonding is detrimental of their con- 
ductive behavior. 



We have performed simulations of the conductance of 
(001) gold and platinum junctions that are bridged by di- 
atomic molecules of the first-row elements, using our code 
SMEAGOL| SMEAGOL is a transport program which 
uses the Non Equilibrium Green's Functions formalism Ji]] to 
compute the charge density and current of the junction. The 
Hamiltonian of the junction is determined by DFT theory via 
the code SIESTaJm]. 

In the case of gold junctions, we have taken for simplic- 
ity a flat surface, oriented the molecules perpendicular to it, 
and set the electrode-molecule distance equal to the equilib- 
rium distance of the gold-H2 junction. We have used a more 
realistic geometry for platinum in contrast, motivated by the 
recent experiments by Tal and coworkers p8|| . We have placed 
in these cases a pyramid of platinum atoms on top of the flat 
surfaces, as in our earlier publications 1 1511 . We have oriented 
the molecules in a variety of angles, including the bridge, per- 
pendicular and tilted orientations. We have finally relaxed the 
forces of the atoms at the pyramids and molecules. We have 
found in this respect that Li 2 , Be 2 , B 2 , C 2 and N 2 relax to the 
bridge position, while 2 , F 2 and CO relax to a tilted orienta- 
tion. 

Our code computes the conductance of the junction via the 
formula 



G(V) 



dIleads{V) 
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d(-eV) 



dio T{E = huj, V) {til — nu) 



where T(E, V) are the energy- and voltage-dependent trans- 
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mission coefficients of the junction, defined in Eq. (17) in 
appendix A, and til,r are the distribution functions of the 
Left and Right electrodes, also defined in Appendix A. G(V) 
can be approximated at low enough voltages by the linear re- 
sponse formula 



(a) 



G(V) ~ G T(E = eV, 0) = G T(E = eV) (2) 



We show in Figs. 1 and 2 the transmission coefficients 
T(E) for all the gold and platinum junctions simulated. No- 
tice that we have taken the Fermi energy of the semi-infinite 
electrodes at equilibrium as the reference energy. Overall, we 
find that the zero-voltage conductance G(0) is of order Go, 
with the only exception of the gold-i 7 ^ junction. Further, all 
the transmission coefficients are rather smooth at positive en- 
ergies, while they show peaks in a range of energies below 
Ep. These peaks do not correspond to molecular states, but 
rather to the d-band conduction channels at the electrodes, and 
are naturally located either a few eV below Ep for gold, or ex- 
tend up to Ep for platinum. We will skip below any further 
reference to this d-band conduction channels, since they are 
not relevant for our purposes. 

A closer look at Fig. 1 shows that T for gold junctions is 
either a plateau (as in Li2, Be2 and F2), or the sum of a plateau 
and a broad resonance (in B2, C2, N2 O2 and C). The low- 
voltage conductances are close to Go for Li2 and F2, while 
they are of order 0.3-0.7 G for Be 2 , C 2 , N 2 and CO, and of 
order 1.5-2 Go for B2 and O2. 

The transmission coefficients of platinum junctions in Fig. 
2 show a slightly more complex structure, but are still pretty 
smooth and show no sign of the narrow resonances that mark 
the tunneling regime. For platinum, L12, and N2 show con- 
ductances of about Go, while Be2 and CO have about half 
a conductance quantum, and B2, C2 and O2 have conduc- 
tances in the range 1.5-2 Go- The conductance of F2 is below 
0.1 Go. Our results for CO agree wit h a p revious simulation 
preformed by Strange and coworkers ll25ll . but fail to repro- 
duce the peak at 1 Go shown in the experimental conductance 
histograms II 1811 . As a side remark, we should point out that we 
have found that the conductance of platinum may increase or 
decrease by even a factor of two depending on the placement 
and orientation of the molecule. 



ANALYTICAL MODELS 

Caroli's model in its simplest form, is depicted schemati- 
cally in Fig. 3. It consists of two identical semi-infinite chains, 
called Left (L) and Right (R) electrodes, which sandwich a 
free-standing atom (M), so that the three of them are initially 
unconnected. The chains have a single orbital per atom of 
energy e; electrons can hop between atoms via the Hamilto- 
nian matrix element t. The three pieces are initially held at 
the same bias voltage, and the common Fermi energy Ep of 
the unbiased chains is taken as the reference energy. Since 
the chains are semi-infinite, they can be regarded as reservoirs 



(b) 



FIG. 3: (a) Two chains made up of P atoms (P — > 00), connected to 
a central atom by a hopping matrix element tju. Each of the atoms 
in the chains has a single orbital with atomic energy e; electrons in 
the chains hop from one atom to the next via the hopping integrals t. 
The central atom has a single orbital of energy em- (b) Two atomic 
chains connected to a diatomic molecule. Each atom in the molecule 
has a single orbital of energy em\ electrons hop between the two 
orbitals via the intra-molecular hopping integral T. 



of electrons, subjected to a given equilibrium chemical poten- 
tial [Il.r, whose energy level population is described by the 
Fermi distribution functions rip^R. The central atom is held 
initially in thermodynamic equilibrium by its contact to a third 
reservoir at the arbitrary chemical potential /jt = eVp- This 
determines the equilibrium distribution of the atom um and 
therefore its population. The whole system is subsequently 
biased by a voltage V so that the chemical potentials of the 
chains is shifted to /j*l,r = ieV/2, while the chemical po- 
tential at the atom can still be left equal to the initial value, 
although it is physically more reasonable to reset it to the av- 
erage between fip and (j,r. Notice that the three pieces stay 
initially in equilibrium since they are unconnected. Later on, 
the central atom is connected to the electrodes by adiabati- 
cally switching on the hopping integrals between them until 
they reach their final value tnj- We assume that the system 
is able to reach an stationary non-equilibrium state long time 
afterwards I22II . where there is a total bias voltage between the 
electrodes equal to V, that induces a finite electron current. It 
is important to stress that the charge population at the atom 
is determined by its contact to three reservoirs, each at a dif- 
ferent chemical potential. Notice furthermore that, while the 
coupling to the chains is controlled by iju, the coupling to 
the third reservoir can not be modulated, so that the flow of 
electrons between it and the atom is completely transparent. 

This model can easily be solved analytically using 
the machinery of the Non-equilibrium Green's functions 
formalism l23l 12411 (see Appendix A for a thorough algebraic 
derivation). Fig. 4 shows the results for the Density of States 
at the atom (DOS) 
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FIG. 4: (Color online) Conductance as a function of voltage (solid 
black curves), and density of states at the atom at zero voltage as a 
function of energy (dashed red curves), for Caroli's model. The nine 
panels displayed correspond to different choices of the parameters 

tM, Sm- 
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FIG. 5: (Color online) Conductance as a function of voltage (solid 
black curves) and density of states at the atom as a function of energy 
(dashed red curves) for the diatomic molecule model. The nine pan- 
els displayed correspond to different choices of the parameters tjvr, 
em- The parameter T is set to 1.5 in all cases. 



where g M M is the retarded Green's function at the central 
atom, which is defined in appendix A. The figure also shows 
the conductance G(V), which has been computed in this case 
performing the numerical derivative of the current through the 
junction, to access the high-voltage behavior. 

Notice that in this model the position of the HOMO/LUMO 
levels with respect to the Fermi energy of the electrodes (taken 
as 0) is given by the atomic energy em, while the hopping 
integral t M denotes the strength of the HOMO/LUMO hy- 
bridization with the electrodes. The hybridization causes two 
important and well-known effects on the bare atomic level 
em- First, it renormalizes it, e. g.: the level changes its en- 
ergy. Second, its broadens it, since now an electron initially 
placed at the atomic level can hop back and forth to the elec- 
trodes, and hence the atomic state acquires a finite lifetime. 
These two effects are readily seen in pm, whereby the initial 
delta-like peak corresponding to the atomic states moves and 
broadens to a resonance. The resonance is sharp if tM is much 
smaller than t, or broad if tM and t are of the same order of 
magnitude. What we wish to stress here is that the extent of 
this broadening also determines the behavior of the conduc- 
tance of the junction. 

Fig. 4 illustrates how the different regimes of a junction 
manifest in the DOS and the conductance and how these 
regimes are controlled by the basic parameters of the model. 
The top three panels correspond to the resonant tunneling 
regime. Notice that this regime is seen even for values of 
the hybridization tM as large as 0.2. This regime is charac- 
terized by a sharp resonance in the DOS, located at the energy 
position of the molecular orbital em- Additionally, sharp reso- 
nances appear in the conductance at different voltages, whose 
maximum height can not exceed Go- The low-voltage con- 



ductance is very small except if the molecular orbital is pinned 
to the Fermi energy, in which case G ~ Go- We stress that 
the molecular level must be fine-tuned to the Fermi energy to 
achieve sizable values of the conductance. It is important to 
stress again that the resonant tunneling regime is achieved by 
values of tM which are smaller, but not much smaller, than t. 
Just a factor of five is enough to place a junction within it. 

The bottom panels in Fig. 4 correspond to highly transpar- 
ent junctions, which corresponds to values of tM — 1. The 
molecular orbital is in this case fully hybridized with the con- 
duction channels at the electrodes. The DOS is very broad, 
and has a width of the order of the bandwidth of the conduc- 
tion channels at the electrodes. More importantly, the conduc- 
tance is very flat and its height is very close to Go- Notice that 
this happens for a wide range of positions of the molecular or- 
bital em- Therefore, if the hybridization between the molec- 
ular orbital and the conduction channels is large, then there is 
no need whatsoever to fine-tune the position of the molecular 
orbital. In other words: when tM — 1> the molecular orbital 
is always tuned, provided it is initially located within the band 
of the conduction channels at the electrodes. This is one of 
the central results in this article. 

The middle panels in Fig. 4 show the situation for junctions 
of intermediate transparency, The DOS reflects this crossover 
behavior, where a resonance has already been developed but 
still has a large width. The conductance is very flat, but we 
point that its height has a stronger dependence on the position 
of the molecular level than was the case of the highly trans- 
parent junction. 



It is fortunate that some Ab initio simulation codes 12 111 al- 
low to make approximate estimates of the physical parame- 
ters that appear in Caroli's model. We have indeed recently 
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performed such estimates for a constriction consisting of plat- 
inum electrodes bridged by H2 molecules. We found that in 
this case t ~ £m ~ 5 eV (e.g.: tj^/t ~ 1), and that the anti- 
bonding state of the molecule was located within the s-d-band 
complex of platinum 1 1511 . We hence ascribed that junction to 
the highly transparent regime, driven by the large value of t]\j- 
Most of the junctions that we have simulated in this article, on 
the other hand, correspond to the crossover regime, where the 
conductance is of order Go, but where its exact value has a 
significant dependence on the exact position of the molecular 
orbital. This is also the case of the Benzene junctions dis- 
cussed by Kiguchi and coworkers! 17 1. 

To make a closer contact with experiments on diatomic 
molecules, we have performed a slight modification of Car- 
oli's model, that we call the diatomic model. In it, the central 
atom is replaced by a diatomic molecule as depicted in Fig. 1 
(b). The nice feature of this model is that on the one hand, it 
accounts for both HOMO and LUMO levels and, on the other, 
it is also easily solved analytically (the algebra is relegated to 
Appendix B). This model can be applied directly to the case of 
the platinum- i?2 constriction, where we showed that the elec- 
tronic conduction is carried by the anti-bonding state of the 
molecule, which is strongly hybridized to the platinum con- 
duction bands, while the bonding state lies slightly below the 
edge of the platinum conduction band and hence does not par- 
ticipate in the chemical bond, showing up as a sharp resonance 
in the DOS ll5ll . As stated above, we estimated that for this 
case ijvf ~ t ~ 5eV (£m/* ~ 1)> while T was slightly larger, 
about 6 eV (T/t ~ 1 - 1.5). 

We display in Fig. 5 the results of the diatomic model for 
the DOS at the left atom in the molecule, and for G(V) for 
the case where T = 1.5, and where the bonding state lies be- 
low the conduction band edge and hence shows up as a sharp 
resonance in the DOS. We again find that when tu is much 
smaller than t, both the DOS and G display resonant behav- 
ior, whereby the conductance is always very small, except for 
specific (usually too large) voltages, or when the anti-bonding 
state is finely positioned at the Fermi energy. On the con- 
trary, when tjvf ~ U the resonance broadens, and the conduc- 
tance has long plateaus with values close to Go- The results 
that compare best to our simulation data for the DOS at the 
molecule in the platinum- H2 constriction (see Fig. 3 in Ref. 
Il5tn correspond to the middle panel in the bottom row. Hence 
we expect that this model can describe correctly the platinum- 
H2 constriction if we use the parameters T ~ 1.5, tM ~ 1 
and em ~ — 1. Notice that the conductance indeed displays a 
value close to Go . 

A final note in this section relates to the large variability in 
the conductance obtained experimentally in BDT junctions. 
We remind that the voltage at which the conductance reso- 
nances depends strongly on the energy position of the molec- 
ular orbital for junctions belonging to the resonant tunneling 
regime. In other words, the value of the measured conduc- 
tance depends strongly on the details of the orbital and its 
bonding for this type of junctions. On the contrary, since 
the conductance has a weak dependence on the details of the 



FIG. 6: (Color online) Geometry of (a) Benzene-dithiolate/gold, and 
(b) Benzene/gold junctions. The drawings show the last layer of gold 
atoms at each electrode and the carbon, hydrogen and sulfur atoms 
of the molecules. Fig. (b) also shows the gold apex atoms attached 
to the gold flat surfaces. 



molecular orbital for highly transparent junctions, we expect 
that the experimental variability must be suppressed for them. 



ANALYSIS OF BENZENE AND BENZENE-DITHIOLATE 
JUNCTIONS 

Benzene and Benzene-dithiolate are conjugated molecules 
whose HOMO and LUMO states correspond to 7r-bonded 
p orbitals, which are delocalized throughout the whole 
molecule. BDT/gold structures are archetypical molecular 
junctions, which are used as templates against which both 
theory and experiments are benchmarked. Experimental data 
show always small conductance values, which range from one 
to three orders of magnitude smaller than Go, depending on 
the experiment. Theoretical results tend to shed larger val- 
ues for G, typically of the order of 0.05 to 0.1 Go. This dis- 
crepancy is usually attributed to a poor description of elec- 
tronic correlations by the theoretical models [26]. BDT/gold 
junctions therefore seem to belong to the resonant tunneling 
regime described in the previous section. Interestingly, Ben- 
zene/platinum junctions are highly conductive lfl7ll . showing 
values of G of the order of Go, which indicates that these 
junctions belong to the transparent regime of Caroli's model. 
It is therefore important to understand better the contrasting 
behavior of these seemingly similar junctions. 

In order to do so, we have simulated BDT junctions at- 
tached to (001) flat gold surfaces. BDT molecules attach via 
the thiol-groups, such that the sulfur atoms are placed in hol- 
low sites at each side, as shown in Fig. 6 (a). We have 
set a sulfur/gold distance of 2 A, such that the junction is 
slightly stretched. We have also simulated Benzene molecules 
attached to (001) gold surfaces. The surfaces are initially 
flat, but we have placed an additional gold atom on top of 
each of them, to provide for a preferred anchoring site. Ben- 
zene molecules bind to these electrodes via the 7r-conjugated 
molecular orbitals, yielding a geometry for the junction such 
as that shown in Fig. 6 (b). We have set a distance of 5.0 A be- 
tween the two apex gold atoms, which is again slightly longer 
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FIG. 7: Transmission coefficients T(E) for the two geometries in 
the previous figure: a) BDT; (b) Benzene.. 

than the equilibrium distance of the junction. We note that 
there is no direct hybridization between gold orbitals at each 
side of the junction for such a distance, which is confirmed by 
negligible transmission coefficients T(E). Upon relaxation 
of forces, we have found that the center of the Benzene plane 
is tilted in such a way to maximize the bonding between the tt 
orbitals and the apex gold atoms. 

The transmission coefficients of both junctions are plotted 
in Figs. 7 (a) and (b). We find that they both show broad 
resonances, which classifies them in a crossover ground be- 
tween the resonant tunneling and intermediate transparency 
regimes. Consequently, the effective hybridization between 
the HOMO/LUMO levels and the conduction channels is nei- 
ther too large nor too small. We therefore expect that the con- 
ductance of these junctions must have a relatively strong de- 
pendence on the details of the junctions. BDT junctions have 
two main peaks peaks, placed 1 eV below and 2 eV above 
the Fermi level, respectively. There is a gap in between, so 
that our simulations yield a zero-voltage conductance of about 
0.1 Go- This height is actually reduced by strong correlation 
effects B26I1 . The placement of the two peaks and the gap is 
due to the fact that the 7r-conjugated orbitals are orthogonal to 
the conduction channels of the electrodes, due to the geometry 
of the junction. The geometry of Benzene/gold junctions on 
the contrary is such that the 7r-conjugated orbitals hybridize 
directly with the conduction channels at the electrodes. Con- 
sequently, the broad resonance is moved, or pinned, to the 
Fermi level. 

The reactivity of the thiol groups is therefore detrimental 
of the conductivity of BDT junctions. In other words, BDT 
molecules attach to the electrodes via the thiol groups, yield- 
ing a junction geometry where the conjugated orbitals do not 
bind to the conduction channels. When the thiol groups are 
missing, Benzene attaches to the electrodes via the less reac- 
tive tt molecular orbitals. While the 7r-gold chemical bond is 
less reactive, it has the virtue of providing a direct hybridiza- 



tion between the HOMO/LUMO orbitals and the conduction 
channels at the electrodes. 



CONCLUSIONS 

We have performed simulations of a number of molecular 
junctions, where (001) gold or platinum electrodes sandwich 
first-row diatomic molecules. We have found that these junc- 
tions are highly conductive, which is manifested both in large 
values of the conductance and in smooth transmission coeffi- 
cients T(E). We have used Caroli's model to argue that this 
is a generic feature of the transparent regime of a junction, 
which is driven by a high hybridization between the delocal- 
ized molecular orbitals in the molecule and the conduction 
levels at the electrodes. 

JF acknowledges conversations with J. van Ruitenbeek, as 
well as his sharing his results with us prior to submission. This 
research has been funded by the Spanish government (project 
FTS2006-12117). 



APPENDIX A: CAROLI'S MODEL 

It is convenient from a mathematical point of view to split 
the total Hamiltonian in two pieces: the unperturbed Hamil- 
tonian of the total system Ho and the perturbation that drives 
the system out of equilibrium Hi. We use the basis of atoms 
in the system, whose 2P + 1 states we denote by \i >. Then 
the Hamiltonians can be written as the following matrices 

( Hl \ / f L \ 

Ho = e M Hi = \ f[ f & . (4) 
V H R J V f R / 

where Hl,r are P x P tridiagonal matrices of the form 

/. . . . ■ • \ 

~ _ . t e±eV/2 t . . 
UL - R - . . t e±eV/2 t . ' { ' 

v • •/ 

and f{ = (...,0,0,t M ) and f R = {t M ,0,0...) are P- 
dimensional vectors. 

The retarded and advanced Green's functions of the unper- 
turbed system can be calculated through the equations: 

[w* -H ]g o H ' A (w)=l (6) 

where cj ± = u ± i5, S being an infinitesimal number. Notice 
that Q are also huge matrices of size (2P + 1) x (2P + 1). 

This large set of coupled equations can actually be reduced 
enormously by gaussian elimination of the atoms in the elec- 
trodes until only three states remain, \L >, \M >, \R >. The 
resulting 3x3 matrix Green's functions are called F to avoid 
confusing them with the conductance, which is denoted by 
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G. Their matrix elements (R/A superindices are henceforth 
dropped when there is no danger of confussion) 



9l 
= | g M 
g R 



(7) 



Here ql.r are the surface Green's function of the electrodes, 
and in our one-dimensional model they are equal to 



R.A 
9lr 



:) 2 



4t 2 



(8) 



To simplify matters, we will take henceforth e = and t as 
the energy unit. We also define a lesser Green's function, that 
carries information of the electron occupation in each piece, 



Then the Keldish formalism |2(A 12311 provides a simple 
recipe to compute the Green's functions of the final system 
in the steady state, 



9ll 9lm 9lr 
9ml 9mm 9mr 
9rl 9rm 9rr 



F< = G R [Gq]~ 



[G 



G J 



9lm 9lr 
9mm 9mr 
9rm 9rr 
(10) 

These equations can easily be solved analytically, with the fol- 
lowing result for the retarded/advanced Green's function, 



( 9 k L 
I 9ml 

\ 9r L 



at 







Fo < = g< 



(9) 



g< 



where each g^ R m( w ) ^ s written in terms of the density of 
states pl,m,r an d the Fermi distribution function ul^m.r of 
each unconnected piece, as gf (u) = 2 ir i pi(ui — eVi) n(u> — 
eVi). 



t M F 



E L (w - e M - 

tM El 

El Efl 



tM El El E^ 

t-M E/j 
tM E_r T,r (co — em — El) 
(11) 



where E l ,_r = t 2 M g L:R and D = uj - e M - El - E fl . Like- 
wise, F < can be written as the sum of the following three 
matrices 



L = 



2-KlpL til 
\D\ 2 



\lj - e M - E fl | 2 
tM — £m — E^) 

Eg ( W - EM ~ E£) 



t M (cj - em - Eg) 



L M 
tM Eg 



t M e4 



Eg) 



R 



2Trip R n R 
\D\ 2 



Ef ( W 



*M ^L 



Eg) t M (^ 



tM^L 
t 2 

£m 



Eg) 



Eg(c 
t M (cj 



- t-M 

£m 



£m 



El) 

' 12 



(12) 



\g M D\ 2 



|ElP 
*m E L 
EgEf 



*m Eg 



t 



M 
yA 



Eg^ 

A 
R 



tM E 



These lesser Green's functions enter the calculation of the 
electronic charge and current. For instance, the charge at the 
atom TYjvf can be written as follows: 




This formula for the charge can be rewritten in terms of the 
conventional expression for the equilibrium state 

N eq = - —lm[g MM ]n L (14) 



plus an explicit non-equilibrium term 




where T L , R = i (Eg R - E£ >fl ) - 2 tt t 2 M p L , R {u> - eV L . R ), 
and likewise Tm = 2nt 2 M pm(^> — eVr). 

The current traversing the left contact is computed via the 
formula 
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which can be written as a sum of the following two contribu- 
tions, 



h 



leads 



latorn.L 



— J du\g MM \ 2 T L T R {n L -n R ) 

Go 

e 



J duj T(E = huj, V) (til — n R ) 

Go f , , R ,2 r L Ta/ 
— / <M5mmI — 12 "A/ 



f 2 



Notice that the current traversing the right contact can be writ- 
ten as the sum of Ii eads plus I a t om ,R where 



t Go f , | R |2 T^Tm 

iatom,R — / UUJ\g MM \ — 72 nM 

e .1 t M 



(18) 



Since these atomic contributions apparently break the conser- 
vation of charge, I a tom.(L.R) are usually dropped (as is the 
second term in Eq. (15)). In other words, the total current / is 
approximated by Ii eads - 



where Hl.r are the P x P tridiagonal matrices described in 
Appendix A, and 



Hi = 



/ T L 0\ 

f I 

f R 

V o o f R oj 



(20) 



After decimating the unwanted degrees of freedom, the unper- 
tubed Green's functions look like 



F, 



R,A,< _ 



R,A,< 






-R,A,< 

9m 






R,A,< 
9 R 



(21) 



where 



APPENDIX B: DIATOMIC MOLECULE MODEL 

The unperturbed Hamiltonian and the perturbation in the 
diatomic model take the form: 





(H L 












Ho = 






£m 
T 


T 

SM 






(19) 




I o 








H R j 


I 



~R,A 
9m 



-T 



-T 

- £m 



(22) 



and .9M = -!/""(<? 



M 



~9m) 



tim- The algebra is more tedious, 
but the final result for the retarded Green's function t 2 M F is 
the following: 



1 

D 



f E L [(w - £m)(w - s M - Er) - T 2 } 
tM El (uj — em — E,r) 
tM T El 

V te l e k 



tM El [lo 
t 2 



- em — El;) 



2 M (lo - em - EhJ 

t 2 T 
1 m 1 

tMTT, R 



tMTT,L TElE_r 
t 2 M T tM TT, R 
^m ( w — £ a/ — El) ^a/ E^ (w — £a/ — El) 
i M E fi (w - e A / - El) E fl [(u> - e M ){u - £m - E L ) - 
(23) 



where 13 = (w - e M )(w - e M - E fl - E fl ) - T 2 + E L E fc 
and the equation for the current traversing the left link can be 
written as 

Ilm =~ J (hjT L T R \G MM , | 2 {n r - n L ) (24) 

where the atomic contribution has been neglected again, and 
Gmm 1 refers to the element (2,3) of the retarded Green's func- 
tion matrix, that connects the two atoms in the molecule. 
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